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We propose an experimentally feasible setup with ultracold alkaline earth atoms to simulate the 
dynamics of U(l) lattice gauge theories in 2+1 dimensions with a Chern-Simons term. To this end 
we consider the ground state properties of spin-5/2 alkaline earth fermions in a honeycomb lattice. 
We use the Gutzwiller projected variational approach in the strongly repulsive regime in the case of 
filling 1/6. The ground state of the system is a chiral spin liquid state with 27r/3 flux per plaquette, 
which spontaneously violates time reversal invariance. We demonstrate that due to the breaking of 
time reversal symmetry the system exhibits quantum Hall effect and chiral edge states. We relate 
the experimentally accessible spin fluctuations to the emerging gauge field dynamics. We discuss 
also properties of the lowest energy competing orders. 
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One of the main motivations of studying ultracold 
atoms in optical lattices is the high extent of experimen- 
tal control. Such systems are very flexible and there- 
fore are good candidates for simulating other quantum 
systems, where experimental control is more cumber- 
some. There is a vast number of proposals where ul- 
tracold quantum gases can serve as simulators of con- 
densed matter, or even high energy physics systems (see 
for instance Ref. [I]). An important example of such 
proposals concerns the recent experimental realization of 
trapping and cooling of ultracold alkaline earth atoms 
which could serve for quantum simulators of high 
symmetry magnetism [5]. Despite of these spectacular 
developments, one of the most important goals of quan- 
tum simulators remains still to be realized, namely the 
simulation of quantum gauge theories, which appear first 
of all in high energy physics, but arise naturally also in 
many areas of condensed matter physics, such as physics 
of frustrated systems, or of high temperature supercon- 
ductors [S]. The main difficulty here is to map the many 
modes of the gauge field to those of an atomic ensem- 
ble. Very recent proposals use mixtures of fermionic and 
bosonic atoms, so that the bosons are the mediators of 
the gauge field [71 1]. 

Here we propose another, somewhat simpler, scheme 
with only a single species of ultracold atoms to simulate a 
2+1 dimensional U(l) lattice gauge theory with a Chern- 
Simons term. Our proposal is based on the observation 
that low energy excitations of certain Mott insulators can 
be described by lattice gauge theories |6j |9j . The Mott 
insulator we consider here is formed by spin-5/2 alka- 
line earth atoms, such as -'^^^Yb, which, as was shown by 
Hermele et al. |10j . can exhibit time reversal symmetry 
breaking, and have a so called chiral spin liquid (CSL) 
ground state in a square lattice. CSL states lack any kind 
of long range order, but due to the violation of time rever- 
sal invariance, they are stable also at low temperatures. 
The fluctuations above the CSL state are described by 



a U(l) gauge theory with a Chern-Simons term arising 
from the chiral (time reversal symmetry breaking) na- 
ture of the ground state [TT] . Here we treat the case of a 
honeycomb lattice; we show that the lowest energy spin 
liquid ansatz with one particle per site that respects the 
underlying lattice symmetries is a CSL. We describe also 
two other lowest energy spin liquid states, that can be 
stabilized in certain situations. We show also how the 
dynamics of the emerging gauge fields is measurable by 
spin correlation functions. Obviously, further motivation 
to study spin liquid phases in a honeycomb lattice with 
ultracold atoms is due to a newly rising interest in related 
challenging problems, like superconductivity in graphene 
[T^HT5] , or other forms of time reversal symmetry break- 
ing that appear for honeycomb and pyrochlore lattices 

[BUT]. 

The main advantage of using alkaline earth atoms is 
that the nuclear spin / decouples from the angular mo- 
mentum of the electrons J, which is zero in the ground 
state, and hence collision processes become spin indepen- 
dent. Therefore spin- J isotopes of alkaline earth atoms 
can be described to a very good accuracy by SU(A^) sym- 
metric model Hamiltonians, where N = 2 / + 1, and we 
take / = 5/2, which is realized e.g. by ^^-^Yb. In an op- 
tical lattice the SU(A^) symmetric Hubbard Hamiltonian 
takes the form 

H=-tY^ (cLcja +H.C.) +\Y1 4aCl/3C»/3Qa, (1) 
{i,j),a i,a.Ji 

where Cia annihilates an atom at site i with spin a G 
{— I . . . |}, t stands for the tunneling amplitude, and U 
for the strength of the on-site interaction. 

In the strongly repulsive regime, U ^ t, the motional 
(charge) degree of freedom of the fermions gets frozen at 
low temperatures leading to a Mott insulator state, and 
the system can be described by an effective SU(6) spin 
Hamiltonian O |9] . The ground state of such a Hamil- 
tonian can be a Neel antiferromagnetic state with long 
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FIG. 1: (Color online) The unit cell for the mean-field ap- 
proximation. 



range order, or a spin liquid state without any kind of 
long range order [18] that can be described using the 
(resonating) valence bond concept [BJ [TWE] . Hermele et 
al. have shown that the Neel order is ruled out by energy 
constraints for SU(A^) Hamiltonians for iV > 3 [lOj . 

Deeply in the Mott insulator region(f7 ^ t) and for 
1/6 filling, there is exactly one particle per site. From 
experimental point of view this is an appealing setup, 
in which undesired 3-body losses can be neglected. In 
this case particle tunneling is forbidden due to the very 
high energy cost of the multiply occupied sites. Only 
virtual hoping is allowed, and the Hamiltonian ([I]) can 
be approximated by 



(2) 



with g ~ —At^/U. This effective Hamiltonian acts in the 
restricted Hilbert space of 1 atom per site; this condi- 
tion is enforced by the local constraints: cJ^Qa = 1- 
The Hamiltonian ^ already exhibits local U(l) gauge 
invariance, since it is invariant under the transformation 
Cia — ^ Cia e*^' , in accordance with the local constraints. 

To study the ground state properties and the low en- 
ergy excitations of the system we decouple the quartic 
terms via a mean-field treatment by introducing the av- 
erage of the pair-correlation defined as 

The effective Hamiltonian takes then the mean-field form 



H^^-9j: E( 



(4) 



Though the mean-field Hamiltonian Eq. Q is already 
quadratic, one still needs to make further assumptions 
about the solution in order to obtain a tractable set of 
equations. We choose a hexagonal unit cell containing 6 
lattice sites, as depicted in Fig. [T] such cell respects the 
original lattice symmetries, and contains as many sites 



as needed to form a SU(6) singlet. As a consequence 9 
independent mean-field amplitudes arise: xi- ■■X9^ and 
the description through the Hamiltonian Q inherently 
describes 6 bands inside the Brillouin zone. With the 
introduction of 6 Lagrange multipliers (one for each site 
of the unit cell) the local constraints can be enforced, and 
then the Mott insulator state with 1 particle per site is 
simply achieved by filling the lowest band. 

By diagonalizing the Hamiltonian Q one obtains the 
ground state | ^P-^'J ) ; using then the self-consistency con- 
dition ([3]) , one arrives at a set of equations for the mean- 
field amplitudes and for the local Lagrange multipliers. 
This set of equations is highly nonlinear and has several 
solutions. Because of the local U(l) gauge structure, so- 
lutions related to each other by the transformation 



Xij ^ Xij ^ 



i{Bi~ej) 



(5) 



are equivalent, and have the same energy and the same 
physical spin wave function, obtained by the Gutzwiller 
projection, i.e., by restricting the solution to the space 
with one particle per site: 



n 



(6) 



The Wilson loop Hj — Ylv Xij calculated for each pla- 
quette Vi is invariant under such transformations, and 
therefore is the same for gauge equivalent solutions. Note 
that due to our construction there are 3 nonequivalent 
plaquettes that can be considered as a 3 sublattice ansatz 
on the dual lattice: the triangular lattice formed by the 
plaquettes. Table |l] shows the three lowest energy mean- 
field solutions. The ground state solutions (first two lines 
of Table |l| correspond to xi = X2 ■ • ■ = Xe = ^ e^'"'^ and 
X7 = Xs = X9 = with r 0.82651 and ip — 27r/9. 
The corresponding Wilson loops are H^ = e^**, with 
$ = 27r/3 flux. Therefore the ground state breaks time 
reversal invariance, and is a CSL state. Due to the vi- 
olation of time reversal symmetry there is a double de- 
generacy in addition to the gauge structure. A higher 
energy staggered like phase is the next-lowest lying quasi- 
plaquette state that has a triple degeneracy with an en- 
ergy of E = —6.062, as it is shown in the following 3 lines 
of Table [ij This phase is the honeycomb analog of the 
(staggered) 7r-flux state of spin- 1/2 fermions on a rectan- 
gular lattice, which is the lowest energy spin-liquid state 
for the non-doped Mott insulator. In our case, due to 
the frustrated nature of the dual lattice, the staggered 
flux phase is not energetically favorable. In this state the 
fluxes do not alternate sign, but they change between 
and TT in a way that every zero flux plaquette is sur- 
rounded by plaquettes with tt fluxes. The last three lines 
of Table |l] show a valence bond crystal type ordering, 
where there is no flux threading through the disconnected 
plaquettes. This state is similar to the ground state of 
e.g. the spin-3/2 fermions in the square lattice at quar- 
ter filling [22] . This valence bond crystal also has a triple 
degeneracy. Figure [2] illustrates these three lowest energy 
states. 
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-6.148 


-0.159 - 0.276j 


-0.159- 0.276i 


-0.159 - 0.276i 


-6.148 


-0.159 + 0.276j 
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-0.223 
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0.460 
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TABLE I: Mean-field solutions. Tlie first column represents 
the energy of the mean-field solution, the other three columns 
give the Wilson loops of the 3 different plaquettes. 
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FIG. 3: (Color onhne) The energy spectrum of the system for 
open boundary conditions in one direction. 



In the CSL state the mean-field generated fluxes pene- 
trating the plaquettes are analogous to those created by 
a magnetic field. As a consequence, quantum Hall ef- 
fect can be observed. It is entirely generated in this case 
spontaneously by the mean-fields. It leads to appearance 
of chiral edge states, i.e. current carrying states localized 
at the boundaries in such a way that the opposite direc- 
tions are separated to opposite edges. Figure [3] depicts 
the energy spectrum with the bulk bands and the edge 
states. For 1/6 filling only the lowest bulk band is filled, 
which is well separated from the next band. We have 
calculated the Chern number (C), which characterizes 
the quantum Hall effect by giving the number of counter 
propagating edge state pairs inside a bulk energy gap. 
We have found that C — 6 for the lowest gap, which 
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FIG. 2: (Color online) Lowest energy mean-field solutions. 
Subfigure a) illustrates the chiral spin-liquid configuration 
with all bonds having the same magnitude. Subfigure b) de- 
picts the quasi plaquette phase with real Wilson loops. The 
dashed line represents a smaller bond value than the solid 
line. Subfigure c) shows the plaquette phase configuration 
with only one nonzero Wilson loop per three cells. Links 
with a zero mean-field value are removed from the figure. 



is in accordance with the observation that there is one 
edge state pair per spin component. It follows that an 
elementary flux of $o = tt/S is attached to every spinon 
(the fractionalized quasi-particle excitation, created by 
the operators cJa)' which in effect renders their statis- 
tics to an anyonic one. On the other hand it provides 
explanation for the <& = 27r/3 flux per plaquette of the 
ground state: since each site belongs to 3 plaquettes, ev- 
ery plaquette contains 2 spinous leading to $ = 2^q, as 
expected. 

The low-energy fluctuations of the mean-field theory 
are given by the phase fluctuations of the mean-field 
amplitude, Xij — Xij^ ^'"'^ , and by the lowest energy 
spinon excitations. The mean-field phase fluctuation 
aij is a gauge fleld with the transformation property 
— > Uij + 0i — 9j. The low energy spinon excitations 
are coming from the vicinity of the energy maxima of the 
valence band, and from the energy minima of the conduc- 
tion band. Since there are 6 such points, the low-energy 
dynamics is governed by spinous with 6 flavors interact- 
ing with a U(l) gauge field. By integrating out the high 
energy (gapped) spinon fields the resulting low energy 
theory is described by the following Lagrangian density 
in the continuum limit: 



An 



ax 



E 

1=1 



(7) 



where the constant arises from the integration of the 
spinon fields and is in the order of the spinon gap g. 
The fields e, and h, are the artificial electric and mag- 
netic fields, respectively, obtained from the scalar and 
vector potentials oq and in the usual way. The con- 
stant C = 6 is the Chern number, and the speed of 
sound V is proportional to 1/g ^ U/t^. In the spinon 
part of the Lagrangian nis stands for the effective mass 
of the spinous, which can be obtained from the curvature 
of the spinon dispersion around the 6 minimas. Due to 
the Chern-Simons term, the gauge bosons mediate only 
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a short range interaction between the spinons, and the 
mean- field solution is stable [9]. Therefore, the low en- 
ergy properties of the system can be characterized by 
U(l) gauge field coupled to the spinons as matter fields. 
The corresponding gauge theory described by the La- 
grangian ([7| can be thus simulated with ultracold spin- 
5/2 fermions loaded into hexagonal lattice. 

Ultracold alkaline earth atoms are produced routinely 
nowadays. A honeycomb optical lattice can also be 
created by sophisticated laser configurations [23|. An- 
other, even cleaner experimental implementation would 
be to use the holographic methods of Greiner et al., or 
Esslinger et al., where an arbitrary two dimensional lat- 
tice potential can be created with the help of an opti- 
cal imaging system [Ml [55]. Detecting the CSL phase, 
or the emerging dynamical gauge theory is not straight- 
forward, but possible. For example, one can measure 
nearest neighbor pair correlations [26], but there is only 
access to \xij\- In fact, according to the Elitzur's theo- 
rem non gauge invariant quantities, such as Xij, average 
to zero [27l [25 . A gauge invariant quantity sensitive to 
chirality and possible to measure is the phase of a loop, 
which can be detected directly by measuring 3-spin corre- 
lations: Si ■ {Sj X Sk)- It's nonzero value witnesses for the 
chiral nature of the spin liquid phase [TT] . Finally, and 
more importantly, one can measure in experiments, using 
for instance spin polarization spectroscopy [29], the spin 
structure factor S^^{i,i]t) — {Sf{t)Sj{0)) at site i and 
J at time t and zero, respectively. This quantity can be 
expressed with the help of the four point spinon Green's 
functions, and in the RPA approximation is given by 



with F^p being the ^-component of the spin matrix in the 
5/2 representation. The quantity n(k,w) is the spinon 
polarization function, which to lowest order is simply the 
contribution of the bubble diagram, and 15(14, is the 
photon propagator up to a numeric matrix for contract- 
ing the space time indices. Due to the hybridization of 
the gauge field and spinon propagators one can find reso- 
nances in the magnetic response function (|8| belonging to 
the spinon and to the gauge field excitations. Therefore 
magnetic response measurements are suitable to reveal 
the gauge structure of low-lying excitations, and the chi- 
ral nature of the ground state. 



We have studied the one particle per site Mott insula- 
tor phases of spin-5/2 ultracold alkaline earth atoms in 
the honeycomb lattice. We have found that the ground 
state is a chiral spin liquid state with broken time rever- 
sal symmetry. Thanks to the finite gap appearing in the 
spinon spectrum, we have integrated out the high energy 
spinon fields and arrived to a dynamical U(l) gauge field 
theory with a Chern-Simons term. This gauge theory de- 
scribes the spin fluctuations of the system, and therefore 
the gauge field dynamics can be accessed experimentally 
with the help of spin response measurements. 
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